Inverse Power Law Quintessence with Non- Tracking Initial Conditions 
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A common property of popular models of quintessence dark energy is the convergence to a common 
solution from a large range of the initial conditions. We re-examine the popular inverse power- 
law model of quintessence (where the common solution is dubbed as the 'tracker') with particular 
attention to the initial conditions for the field and their influence on the evolution. We find that 
previously derived limits on the parameters of the potential in this model are valid only in a range 
of initial conditions. A reasonably sharp boundary lies where the initial energy density of the scalar 
field is equal to that of the background radiation component. An initial quintessence energy density 
above this equipartition value lead to a solution that will not have joined the tracker solution by the 
present epoch. These non-tracker solutions possess the property that their present equation of state 
is very compatible with the observed bounds and independent of the exponent of the potential. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 



Modern observations favor a dark component of 
the cosmological energy density responsible for the 
present accelerating expansion rate of the Universe 0. 
Quintessence is a proposed model for this component 
in which a scalar field slowly evolves towards the min- 
imum of its potential @ [j As the quintessence field, 
Q, nears this minimum and the density of matter, pm, 
simultaneously decreases, the expansion rate will natu- 
rally start to accelerate. Making Q dynamic alleviates 
the large discrepancy between the energy densities of ra- 
diation, matter and dark energy, most evident during 
the early Universe. The evolution of such a scalar field 
depends entirely on its potential V{Q). Several poten- 
tials have been proposed and investigated, two in partic- 
ular have gained interest because of their simplicity: the 
inverse-power law (IPL), V{Q) = U Q" 7 llJilllll3, 
and the exponential, V(Q) = e~ XQ , [H Illf7 Both of 
these potentials possess attractor-like solutions such that 
a broad range of initial conditions converge to a unique 
solution at late times 0, 0, . Typically, during its 
evolution the ratio of the energy density in Q to the criti- 
cal density of the universe, Oq, at a fixed redshift is larger 
than the energy density of the pure cosmological constant 
it hopes to replace, regardless of the potential used. In 
some cases, for example the exponential potential, the 
energy density during the early Universe may be large 
to affect Big Bang Nucleosynthesis (BBN) [Tol ITil lit). 
But in other models, such as the IPL, pq is a small and 
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slowly increasing contributor of the energy density until 
close to the present time when pq ~ pm- 

Probes like BBN, the Cosmic-Microwave Background 
(CMB), Type la Supernovae (SN la), and Large Scale 
Structure (LSS) place restrictions on the parameters ap- 
pearing in a potential. By combining CMB, SN la, and 
LSS data, Bean & Melchiorri obtained a limit on the cur- 
rent equation of state of the field, wq, of —0.85 > wq > 
— 1 (at la) ^5|. In the IPL model these results led to 
the conclusion that the exponent for the power law, 7, 
must be < 2 at 95% confidence 16] . Similar results were 
found in Corasaniti and Copeland |17j , leading these au- 
thors to all but rule out IPL models. Doran, Lilley, and 
Wetterich 18] combined CMB and LSS data and found 
that IPL models with 7 < 1 were favored while Yahiro et 
ai ^3 determined that BBN can constrain tracker mod- 
els with a conservative bound on the equation of state, 
wq < —0.2, and a more stringent bound on the initial 
energy density, pq/pb < 0.01. Yahiro et al note that 
BBN provides a better constraint for models in which Q 
does not track at this epoch, limiting larger pq in these 
models. 

All these previous limits are derived from the assump- 
tion that Q has reached its tracker solution before the 
cosmological test is applied. In this paper we examine 
the range of initial conditions of the IPL potential that 
lead to an evolution for Q that does not converge to the 
tracker solution by the present epoch. We begin by dis- 
cussing the IPL model for Quintessence and its tracker 
solution in section [D] with emphasis on the initial condi- 
tions of this solution. In section flllll we proceed to exam- 
ine non-tracker solutions and the behavior of Q when the 
initial conditions vary significantly from the tracker so- 
lution. Finally, in section Hvl we re-derive the constraints 
on the parameters in the potential by using the CMB 
and SN la observations when the tracker initial condi- 
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tions are used, then we will examine how the constraints 
on the potential parameters are altered when the initial 
conditions are altered. 

Throughout this paper we will assume a flat geome- 
try for the Universe based on the location of the first 
peak in the CMB angular power spectrum [2£j, EI and 
work in reduced Planck units, Mpl = 1/vo7tG = H = 
c = 1. Whenever we require numerical values for the 
cosmological parameters we choose Qm — 0.3 and 
H = 72km/s/Mpc 



II. TRACKER QUINTESSENCE 

Our focus in this paper is initial conditions and how 
they affect of Quintessence at late times. Since we will 
frequently make comparisons of the solution calculated 
from arbitrary initial conditions with the 'tracker' solu- 
tions we begin with a derivation of the properties of the 
latter. 

As with all the non-interacting fluid components in the 
Universe, the rate of change of the quintessence energy 
density, pQ, obeys 



through J5J two first order, coupled, equations 



PQ = -3H(p Q + P Q ) 



(1) 



where Pq is the pressure of Q, the dot represents differ- 
entiation with respect to cosmic time t, H is the Hubble 
parameter, 



a at 



PQ + Pf 



(2) 



a is the scale factor and pf is energy density of all the 
other fluid components. From this equation, and the 
definition of the energy density and pressure as 



PQ = \ + V{Q) 
Q 2 

Pq = \~V{Q), 



the equation of motion for Q is simply 
Q + 3HQ + V = 0, 



(3) 
(4) 

(5) 



where V is the derivative of the potential with respect to 
Q. From these general equations there are many different 
paths that would allow us to advance, each emphasizing 
different properties of the field, the solution or the po- 
tentials (see, for example, 0, Q, El). With a general 
potential they are all equivalent of course. For us the 
direction we take is to introduce the equation of state 
wq = Pq/ Pq as the ratio of the quintessence pressure 
to its energy density and a = Pq/ Pf as the ratio of the 
quintessence and remaining fluid energy densities. For 
these two quantities we can derive from equations 



da 
da 
dwQ 
da 



3a(wf - Wq) 
a 

3(l + w Q )(w Q -v 2 Q ) 



(6) 
(7) 



where Wf = Pf / Pf is the equation of state for the re- 
maining fluid components, a is the scale factor, and Vq is 

simply Vq — Pq / pq . The density parameter Dq , the ra- 
tio of the quintessence total energy density to the critical 
density, is related to a in our flat Universe by 



Q 



1 



(8) 



These equations are valid for any potential V(Q). 

For radiation and non-relativistic matter, w and v 2 are 
constants, but for the field this is not the case, and some 
remarkable properties of quintessence can be attributed 
to this fact. If we now adopt the inverse power law form 
for V 



V=*L 
Q 7 



(9) 



then the speed v 2 q in equation J7J is the rather compli- 
cated 




3 Q H 2 {l-w Q ) 



2V 



l/ 7 



(10) 

Obtaining a solution for the evolution of the field is then 
a simple case of integrating equations © and J7J, with 
specified values for 7 and Vb and initial values for a and 
wq. These four parameters are not, however, indepen- 
dent because the present state of the Universe is implicit 
in H and Wf. To calculate H we need to specify the en- 
ergy densities of the fluid components at all epochs: for 
the quintessence we can simply construct this quantity 
according to equation Q, while for radiation and mat- 
ter this can be done by their simple scaling with a along 
with knowledge of their present energy densities. The 
contributions to the current radiation density include the 
measured temperature of the CMB, and the calculated 
energy density for three, non-degenerate neutrino flavors. 
This leaves only the present matter density, which is usu- 
ally expressed as the product of its density parameter m 
and the present critical density 3Hq. If we impose the 
requirement of flatness then we must have Qq = 1 — CIm 
at the present time. We have thus specified the value 
that £Iq must reach at the end of our integration thereby 
introducing a constraint and eliminating one degree of 
freedom. If we regard the initial values for a and wq as 
inputs, and take 7 (V"o) as the remaining free variable, 
then we must adjust Vq (7) to produce a self-consistent 
cosmology. 

The equations 10 and (J7J) for a and wq are clearly 
non-linear but certain cases are open to analysis. In par- 
ticular, we can look at the situations when either deriva- 
tive is zero and then examine whether these conditions 
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are stationary. From equation JBJ we find that trivially 
da I da — when a = and obviously this is stationary 
since all higher derivatives of a with respect to the scale 
factor vanish. There is also a second solution that oc- 
curs when vjf — wq independent of the value of a. If Wf 
is a constant then to be stationary we require wq to be 
also independent of the scale factor but, as we will show 
below, this only occurs in the limit 7 — > 00. 

For equation Q, dwQ/da = occurs when wq = ±1; 
the wq = — 1 case is clear from equation {7J|, while the 
wq = +1 solution arises because Vq = +1 in such cir- 
cumstances. These extrema for the equation of state are 
only stationary for certain potentials because they re- 
quire either constant zero kinetic energy (wq = —1) or 
potential energy (wq = +1). For the IPL they are impos- 
sible to achieve exactly but wq = ±1 can be approached 
in the limit when pq is dominated by either Q 2 /2 or V. 
There is also a third solution that occurs when wq — Vq , 
and neither has unity magnitude. The equality only oc- 
curs for specific values of a and wq and to be stationary 
this third solution requires a balance between the scale 
factor dependence of the \ZTIq and [Qq H 2 } 1 / 1 terms in 
equation l(TUjl. i.e. 



= constant. 



(11) 



This result is equivalent to equation (5) of 6]. Inserting 
the definition for a this transforms into 



a 



1+2/7 



1 



1 + ai 



6 (l+wf)/f 



(12) 



with the subscript / denoting initial values. But if wq 
and Wf are constant then from equation © we also have 



ai 



ai 



3 {w;—wq) 



(13) 



Thus we have two equations involving a, (|12fl and (|13fl . 
which must be self consistent in their scaling with a. 
Taking the limit 7 — > 00 in equation l|12l) leads to a 
varying very slowly with the scale factor, and therefore 
from equation (|13fl . wq —> Wf. This result validates our 
previous statement regarding the stability of da /da = 0. 
Alternatively, taking the limit a, ai <C 1 in equation l|12|) 
and equating it with equation (|13fl produces the require- 
ment that 



wq 



JWf 



7 + 2 



(14) 



for this case. As noted by Steinhardt, Wang and Zlatev 
0- wq < Wf so that from equation (|13f) we see that a 
increases with the evolution of the Universe. These solu- 
tions where dw/da = and a <C 1 are the 'tracker' solu- 
tions originally named by Steinhardt, Wang and Zlatev 
and Zlatev, Wang and Steinhardt 0. As shown in 
IT^ . these solutions are stable attractors for all posi- 
tive values of 7 in that perturbations away from this so- 
lution are damped. As we stated explicitly, this epoch of 




constant wq only occurs when Wf is also constant so if we 
allow more than one component for the background fluid 
then whenever the Universe transits from domination by 
one component to another there will be a period during 
which dwQ/da ^ 0. Finally, if wq is constant during an 
epoch, the energy density of Q scales as pq cx a~ 3 
from which, after using the definition of the equation of 
state, equation l|14fl and the scale factor-time relation, 
we can obtain 



(> = (>! 



2/(7+2) 



(15) 



where Qi is the value of the field at the initial time ti. 
This is the result obtained by Liddle and Scherrer |13| . 

For quintessence to be a viable explanation for the cur- 
rent epoch of acceleration, it must have recently come 
to dominate the energy density of the Universe, i.e. the 
Universe is transiting between matter domination to dark 
energy domination. In such circumstances the analytic 
solution in equation (|15fl is no longer valid and the equa- 
tion of state is no longer a constant. But there is still a 
solution to which all others converge and we shall con- 
tinue to refer to this solution as the 'tracker' even if the 
evolution of Q is no longer given by equation (|15f) and the 
name no longer justified in its original sense. The appeal 
of IPL models of quintessence is not that the analytic 
equation for the tracker is always valid or that dwQ /da 
is always zero but rather that a wide range of initial con- 
ditions converge to a common solution thereby avoiding 
any fine tuning of the initial values of a and wq. We shall 
discuss further how the tracker is approached in the next 
section, for now we will assume that the initial conditions 
are exactly those of the tracker solution. 

The tracker solutions correspond to specific initial val- 
ues of a and wq, the latter stems from equation i|14|) 
while the former can be derived from the requirement 
wq = Vq and equation 11U|) . Since the initial conditions 
are fixed the tracker solution is a function of only one 
parameter, 7. We start integrating the equations at an 
initial scale factor a/ of aj — 10~ 30 : this initial value 
for the scale factor is, of course, arbitrary and we have 
checked that our results are not sensitive to its value. 

Figure shows the initial tracker value of a at our 
choice for the initial scale of aj = 10~ 30 as a function of 7 
The figure shows that ai grows rapidly as 7 increases so 
that even modest values of 7 raise the initial dark energy 
density by many orders of magnitude above the cosmo- 
logical constant value. In figure © we have plotted the 
equation of state at the present time as a function of 7, 
together with the value during both matter and radia- 
tion domination as given by equation (|14fl . This figure 
can be used to derive an upper limit of 7 based simply 
on the requirement that the Universe be accelerating at 
the present time. To accelerate in a flat Universe we 
must satisfy wq Qq < —1/3 so with our adopted value of 
£Iq = 0.7 we find 7 < 5. If we use the Bean & Melchiorri 
3 — (7 limit of wq < —0.55 (though this was obtained by 
assuming wq is constant) then the limit drops to 7 < 3. 
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FIG. 1: The value of the ratio of the quintessence and back- 
ground fluid energy densities, ai, at the initial scale factor 
ai = 10~ 30 as a function of the power law index 7. 




FIG. 2: The equation of state, wq, at the present time (solid), 
during matter domination (long dashed) and during radiation 
domination (short dashed) as a function of the power law 
index 7. The three horizontal, dashed lines correspond to the 
1, 2 and 3-cr limits of Bean & Melchiorri 151 



The equation of state of the tracker at the present time 
is somewhat sensitive to the adopted energy density f^Q . 
Larger Qq leads to slightly smaller wq for a given expo- 
nent so that enforcing wq £Iq < — 1/3 leads to a limit on 
7 that rises slightly more quickly with £Iq than figure © 
would seem to indicate. 



III. NON-TRACKING QUINTESSENCE 



the focus of this paper, it seems reasonable to ask what 
are the 'natural' initial conditions for Q. 

Malquarti & Liddle followed a quintessence field 
with an IPL potential during inflation and determined 
that due to quantum fluctuations Q will diffuse to a very 
low energy density compared to the background fluid at 
the end of inflation. In this scenario Q is decoupled from 
all other fluid components, including the inflaton. Dur- 
ing inflation, the small scale set by Vq in the Q potential 
implies that the inflaton field dominated the dynamics of 
the universe. Only for small values of Q could the dynam- 
ics of the universe be affected by quintessence. However, 
these authors found that at the end of inflation, pq is 
initially so small that the field does not reach the tracker 
until the relatively late stages of cosmic evolution. 

Alternatively we have initial conditions for Q moti- 
vated from inflation where the scalar inflaton field, <fi, 
transits through a period of kinetic energy domination 
after inflation. The epoch in which the kinetic energy of 
the scalar field dominates the expansion rate has been 
called 'kination' [24|]. Previous authors have attempted 
to connect the scalar field in quintessence to the field 
responsible for inflation, by transiting the field through 
a kinetic energy dominated regime |2^, [2(| [53, H^- If 
the scalar field, here responsible for both inflation and 
quintessence (i.e. cf> = Q), does go through a period of ki- 
nation, this period must end by BBN so that the universe 
is radiation dominated at this stage (2t|. This allows us 
to constrain the energy in Q after inflation. Starting 
from the scaling relation for temperature, T ~ a -1 , as 
well as the scaling relations p r ~ a -4 and p^ ~ a~ 6 
during kination, and expressing the relativistic energy as 
p r = ir 2 g*T /30, where counts the total number of 
relativistic degrees of freedom, we obtain 



30 T\ 2 



(16) 



for the energy in <j> at the end of inflation. Here Ti 
and T* are the temperatures at the end of inflation 
and at the end of kination, respectively, and 5* ~ 100. 
As noted in |3(j, in order to satisfy the nucleosynthe- 
sis bound, the temperature of particles created during 
kination just after thermal equilibrium is achieved is 
T eq > 10 6 GeV. Since T eq < T u we can take 10 6 GeV 
to be a very generous lower bound for Tj, which in l|l(j|) 
gives p^.i > 10 42 GeV 4 . Thus it seems prudent to con- 
sider initial conditions on Q such that pq ~ p r in the 
early universe. 



Initial Conditions for the Q-field 



Approaching the Tracker 



Most previous studies of cosmological constraints on 
the IPL model potential have assumed that Q reached 
the tracker long before the epoch at which the tests are 
applied. This is justified since the appeal of the IPL 
model is the fact that Q reaches this solution for many 
initial conditions. Since the initial conditions for Q are 



In has long been known that the tracker solution in 
IPL models (by which we mean the solution given by 
equation (|15fl when flQ is small or what this solution be- 
comes when J1q is large) is an attractor but that under 
certain initial conditions the time required to reach the 
tracker solution will be longer than the present age of 
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the universe |6j, |31| . Several authors have derived limits 
on 7 based on the argument that the tracker must be 
reached at some point before the present epoch. For ex- 
ample: Steinhardt, Wang & Zlatev [(| bounded 7 > 5 by 
demanding that the field reaches its tracker solution by 
matter-radiation equality when initially the energy den- 
sity of the Universe at the inflationary scale was equipar- 
titioned between the Quintessence and the remaining 
fluid components. Similarly, Frieman and Waga Q ob- 
tained the bound 7 > 3.6 using SN la data for currently 
acceptable values of the matter density, fijvf ~ 0.2 — 0.4. 
We shall not impose these restrictions and consider in this 
section all solutions that are able to satisfy the boundary 
condition whether they reach the tracker by the present 
epoch or not. Since we shall frequently compare a so- 
lution calculated from arbitrary initial conditions to the 
tracker we will denote any quantity associated with the 
latter by the subscript T. 

We introduce the quantity u, as in defined to be 
the ratio of the solution Q{a) given arbitrary initial con- 
ditions to the tracker solution Qx(o-) for some fixed value 
of the IPL exponent 7 i.e. 



u(a) 



Q(a) 

Q T ( a y 



(17) 



With this definition of u(a) a solution reaches the tracker 
solution when u = 1, du/da = 0. When fix <C 1 we can 
derive an equation of motion |l2j | from equation JSJ after 
using equation (|T5l) 



rn + 



7 + 2 Wf + l 



7 



t ii 



7 + 2 



7 + 2 w f + 1 



(18) 
= 0. 



In figure J3J) we plot u(a) for the case 7 = 4 for four 
values of ai, the value of a at the beginning of the in- 
tegration, from a = 1, 10~ 25 , 10" 75 , lO" 100 . In all cases 
we start the field from rest and at an initial scale factor 
of a 1 = 10~ 30 . From figure we see that the tracker 
solution for 7 = 4 corresponds to a value of ax ~ 10~ 40 
at ai = 10 -30 so our demonstration values span a large 
range around the tracker. Despite this large range for ai, 
both above and below ax, the figure shows that the solu- 
tions for all are remarkably similar in their general shape 
and in every case we have the remarkable result that the 
tracker solution is always approached from above. The 
oscillations around Qx seen in Liddle & Scherrer are 
too small to be seen with the scale of the figure because 
they occur when u ~ 1. For the two cases where ai < ax 
the initial value of Q is much larger than the tracker Qx 
thus the initial u is much greater than unity. As the Uni- 
verse evolves u falls as a power law with a and the two 
solutions join the tracker at approximately a scale factor 



of a - 10~ lb for on = 10 



-75 



10" 7 for a T = 10- 



100 



In contrast, when ai > ax the initial value of u is smaller 
than unity. As we integrate the equations we see a very 
rapid and immediate increase as u rushes up to a peak 



value after which u turns over and again falls as a power 
law. For the two cases here the solution joins the tracker 



at a 



10- 



for ai = 10 



1 for a 1 = 1. From 



these four cases it is apparent that a value of ai in rela- 
tive proximity to the tracker solution ax joins the tracker 
very quickly but as ai strays from ax the time required 
to reach the tracker grows rapidly. The figure clearly 
shows that the epoch at which the tracker is joined can 
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FIG. 3: The ratio of the value of Q to the tracker as a func- 
tion of scale factor. The four curves correspond to different 
initial values of ai, the ratio of the initial Quintessence energy 
density to the background (radiation) component. In the top 
panel, from top to bottom, ai € {1, 10 -25 }, while in the lower 
panel, from top to bottom, en G {1CT 100 , 1CT 75 }. 
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FIG. 4: The evolution of wq as a function of scale factor 
for the values of ai in figure In the top panel, oti = 1 
(solid) and oti — 10 -25 (dotted) while in the lower panel. 
ai = 10" 100 (dotted), ai = W' 75 (solid). For ai = 1, the 
tracker is not reached by the present epoch and in this case 
wq is only beginning its climb to the tracker presently. 



easily be after the epoch of BBN, of recombination, or 
indeed it may still not have happened. 

In figure (@J we plot wq versus scale factor for 7 = 4 
and the same four values of ai . The behavior of wq in 
non-tracker phases for Q within the IPL model was also 
discussed by Zlatev, Wang and Steinhardt 0. This fig- 
ure helps explain the results in figure 10) • Every case 
is initially set at wq = — 1 and for aj < cyt the equa- 



tion of state remains at this value until the tracker is 
joined at a ~ 1CP 16 and a ~ 10~ 7 seen previously. Note 
that in this figure the oscillations of the field around the 
tracker are seen as the small wiggles in wq immediately 
prior to the joining of the tracker solution. During this 
phase of wq = — 1 the field behaves as a cosmological 
constant. The figure indicates that the larger the dis- 
crepancy between aj and ar for aj < ay the longer 
this pseudo-cosmological constant phase and eventually 
we would expect that at some sufficiently small a] the 
solution never joins the tracker. However, such a small 
value of aj circumvents the main appeal of Quintessence 
and we will not pursue this possibility further. In con- 
trast, for those cases where a > ax the equation of state 
jumps from wq = — 1 to wq = +1 very rapidly and re- 
mains there for some time before promptly plummeting 
to wq — — 1. This period of kinetic energy domination 
is much longer than the rise and turnover seen in fig- 
ure © and its duration grows as the difference between 
0.1 and ar increases. Once again, during the phase of 
wq = — 1 the field behaves as a cosmological constant 
and from there the solution joins the tracker in the same 
fashion as the two cases where a < olt- Eventually we 
expect that at some sufficiently large ai the equation of 
state never falls to wq = — 1 and so the limits on wq ex- 
clude values of aj larger than this. From our numerical 
simulations we find that wq > at the present time for 
a i ^ 10 83 when 7 = 4. 

The solution for u in the limit u 3> 1 is very different 
from that of Liddle & Scherrer [12|. When u ^> 1 the 
l/it 7+1 term in equation (|18|) is negligibly small and can 
be dropped. The solution of the equation is then the sum 
of two power laws 



t 

U 

where the exponents are 

Wf-1 2 



+ a- 



A± = 



2 («;/ + !) 7 + 2 



- 1 



2(117 + 1) 



(19) 



(20) 



Since Wf < 1 the two exponents are both negative and 
A+ > A_. In terms of the initial conditions the two 
constants, u+ and u- , are given by 



ui 



Wf 



1 



Wf - 1 
Wf + 1 



UI ti 



2ui 
7 + 2 



Wf 



1 



Qi t i 
Qi 



■ UI 



(21) 



(22) 



From the second result in equation \\1'2\ we see it_ < 0. 
At late times only the first term in (|19fl is important 
and, after referring to equation H15|) , we see that its time 
dependence is that of the tracker. The field Q is fixed. If 
we start the field from rest then u_ = and so u + = uj. 
This is the behavior see in the bottom panel of figure 10 • 
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From examining the top panel of figure J2J we see that 
empirically this same formula also applies to the cases 
when a i > olt- Thus we are led to conclude that Q is 
essentially stationary even though the field may be dom- 
inated by its kinetic energy. Such behavior is expected 
since nothing in the derivation of equation Ijl9|l depended 
upon wq, the only requirement was that u 1. How- 
ever, in these cases the initial value of u is much smaller 
than unity so strictly equation ljl§|) does not apply ex- 
actly at ti and u + and u_ are therefore not given by 
equations l|21[) and 1)22(1 . But this limit is reached very 
quickly and, furthermore, the peak value of u is approxi- 
mately the same as u+, the normalization of the late time 
evolution. By locating the time t* where the derivative 
of u vanishes we find that 

^t^—j—. (23) 

where the subscript * indicates all quantities are eval- 
uated at t*. First, we can obtain an estimate for Q* 
by assuming the Quintessence energy density at the end 
of the initial ascent is approximately equal to the ini- 
tial Quintessence energy, pi, and, furthermore, that it 
is dominated by its kinetic component. Under these ap- 
proximations Q* « y/2 pi. Second, the peak is located 
so close to the beginning of the integration that £* w ij 
and Qt* ~ Qti so if we make use of tj ~ 1 /Hj and the 
definition of H the result is that 

where f2/ is the initial Quintessence density parameter 
and we have introduced the quantity Q*. 

Before proceeding to further investigate the implica- 
tions of the non-tracker solutions we temporarily digress 
to examine the importance of the initial scale factor upon 
our results. So far we have begun all our integrations of 
a and wq from an initial scale factor of ai = 10 -30 . 
We already know the tracker value of a increases with 
the scale factor from equation (|13|l and as we showed in 
figures © and (QJ the closer the value of aj to olt the 
shorter the non-tracker phase. From these considerations 
it is therefore natural to expect solutions with an arbi- 
trary value of ai to reach the tracker more rapidly as ai 
increases. Figure (J5J, a plot of u versus a when a/ = 1 
and for different oj, shows that this is exactly the case. 
As ai increases the duration of the non-tracker phase 
and the height of the peak value of u decrease. But the 
figure also shows that the scale factor at which a solu- 
tion reaches the tracker is independent of a/. This is not 
totally unexpected: our estimate of ~ \J &>i doesn't 
change with ai if we use the same value of ai so the 
change in it* as ai increases simply reflects the evolution 
of the tracker solution. From this result we see that the 
evolution of the field as we approach the current epoch 
is robust and the necessity of specifying the initial scale 
factor for non-tracker solutions has disappeared. 
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FIG. 5: The evolution of u as a function of scale fac- 
tor for 7 = 4, ai = 1 and initial scale factors of a/ 6 

{io- 30 ,io- 25 ,io- 20 ,io- 15 } . 

When 7 = 4 the two cases ai — 10~ 25 and ai = lO -75 
are sufficiently close to ax that the tracker is reached dur- 
ing the radiation dominated era, while the a] — 1CP 100 
curve reaches the tracker at the beginning of the tran- 
sition to the matter dominated era. The ai — 1 curve 
does not reach the tracker at all, and were it not for 
the initial climb to wq = +1, this case would be iden- 
tical to a cosmological constant. This is consistent with 
our expectations if were approximately the same as 
the present value of Qt- From Malquarti & Liddle 0] 
Vq ~ 1CP 120 (in reduced Plank mass units) for 7 = 4 
so together with wt ~ — 1 for the current equation of 
state of the tracker and our concordance values of £Iq 
and Ho then Qt ~ C(l)- Thus, Q* must also be of 
this order implying that aj > 1 in order for the solu- 
tion to avoid joining the tracker by the present epoch. 
Identifying aj ~ 1 as dividing the cases where the tracker 
was joined from those where it has not is confirmed (and 
is much more obvious) in figure 10 where we plot the 
equation of state at the current epoch as a function of 
a 1 for different values of 7. This figure clearly shows the 
change at aj ~ 1 where the equation of state of Q drops 
to wq = — 1 as ai passes through the equipartition value. 
There is a slight movement of the transition point with 
7 but it is not large. 

In figure J7J) we summarize our discussion of non- 
tracker Quintessence with a contour plot of the final value 
oiwQ as a function of 7 and a/. For values of Or < 1, the 
final value of wq is only a function of 7 while for ai > 1 
we have the situation where the equation of state at the 
present time is close to —1 for a large range in 7. The 
slight dip in the contours of wq — —0.5, —0.6, —0.7 are a 
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FIG. 6: The final value of the equation of state as a function 
of cti. The slight jump at qj = 1 is a result of the field 
oscillating about the tracker today. 
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FIG. 7: Curves of constant wq at the present epoch as a func- 
tion of 7 and ai. From top to bottom, the curves represent 
wq = -0.5, -0.6, -0.7, -0.8, -0.9. 



result of Q just reaching the tracker at the present epoch. 
Just before Q 'hooks' the tracker it oscillates around that 
solution and these oscillations are seen as dips in the 
contours. The figure shows essentially a dichotomy of ac- 
ceptable solutions in the IPL model. In order to obtain 
an equation of state within the current observationally 



favored limits —0.85 > wq > -1 we are forced to small 
values of 7 if aj < 1 but if the initial conditions are such 
that the quintessence energy density is comparable to, or 
larger than, the radiation component virtually any value 
of 7 is permitted. A Universe initially dominated by a 
quintessence field evolves to one that effectively possesses 
a cosmological constant, independent of the value of 7. 



IV. QUINTESSENCE, THE CMB, AND SN IA 

We have already derived an upper limit to 7 just from 
the requirement that the Universe accelerate. We can im- 
prove upon this result by using observations of the CMB 
and SN la, which provide powerful probes of cosmology 
from a redshift of ~ 1000 and lower. 

The significant changes introduced into the CMB an- 
gular power spectrum by replacing a vacuum energy with 
an IPL Quintessence field have been discussed by many 
authors including Ratra & Peebles 3 and Kneller & 
Steigman 14]. The position of the peaks in the angu- 
lar power spectrum are proportional to the ratio of the 
sound horizon at recombination to the angular diameter 
distance at that redshift, while the height of the peaks is 
set by their variance relative to the variance at the largest 
angular scales as measured by COBE. For a given redshift 
and matter density, Quintessence increases the dark en- 
ergy density, leading to a swifter expansion rate. For the 
IPL tracker solution, the contribution from pq is small 
until the recent epoch, thus its influence only becomes 
significant at low redshifts. The swifter expansion results 
in a reduction of the angular diameter distance, 

d A {z) = (1 + z)- 1 f , dz'/H{z') (25) 
Jo 

while leaving the sound horizon at recombination unal- 
tered. An increase in H(z) therefore moves the peak 
positions to larger angular scales and they become more 
tightly bunched. This effect increases with 7 because for 
tracker solutions the energy density is correlated with 
the exponent: smaller 7 corresponding to smaller pq 
at any fixed redshift. At the same time models with 
Quintessence always pass through the matter-dark en- 
ergy transition at a higher redshift even if they typically 
begin to accelerate at a lower redshift. The earlier domi- 
nance of the dark energy will result in a more significant 
decay of the gravitational potentials so that any pho- 
ton traversing a potential well will emerge with a larger 
blueshift relative to the background expansion. The in- 
creased blueshifting (and the increased redshifting from 
potential hills) leads to an enhancement in the variance 
of the CMB temperature. This integrated Sachs Wolfe 
(ISW) effect is most effective at scales above the causal 
horizon at the time of the matter-dark energy transition 
i.e. the largest angular scales. But the COBE normal- 
ization of the CMB power spectrum is exactly at these 
scales so the net result is that it is the variance at the 
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FIG. 8: The x P er degree of freedom for the SN la (dotted) 
and CMB (solid) data, with Q in the tracker solution. The 
'goodness-of-fit' limit for the SN la data at 99% is xt = 1-50 
while for the CMB data the 99% limit is at xt = 1-60 

smallest scales that is reduced. Once again, this effect 
increases with 7 for exactly the same reason. 

The recent evidence for the current acceleration of the 
universe comes from the use of SN la as 'standard can- 
dles' at high-redshifts, z ~ 0.3 — 0.8 Q. SN la are good 
standard candles because of the strong correlation be- 
tween the peak in the absolute magnitude and the de- 
cline in brightness 15 days after the peak [3^, IH, l34| . 
A measurement of the apparent magnitude is all that 
is needed to determine the luminosity distance di. This 
same quantity can be computed for a given cosmology 
through 

*-< 1 + *>jfi£r (26) 

From this equation it is immediately obvious that an in- 
crease in the Hubble parameter results in a smaller lu- 
minosity distance and hence a brighter SN la for a given 
redshift. For tracker solutions, as 7 increases the lumi- 
nosity distance decreases and so the same SN la will ap- 
pear brighter still. 

In figure (JBJ we compare the IPL model predictions 
using tracker initial conditions with the CMB and SN la 
data by plotting xt: the x 2 P er degree of freedom v as 
a function of 7. For the CMB data we use the 41 points 
of BOOMERANG jH, DASI ^ and MAXIMA ^ 
together with RADPACK [H] to calculate the x 2 - We 
optimize over the scalar tilt and the baryon-to-photon 
ratio but, in order to make the calculation computation- 
ally feasible, we have not marginalized over the present 
Hubble parameter or the matter density Ojvf ■ For SN la, 



FIG. 9: Contours of constant xt m the 7 — a plane. From 
top to bottom, the curves are: xt ~ 1-60 (solid) and 1.40 
(long dashed) corresponding to a probability of 99% and 95% 
respectively for the fit between theory and the CMB data, 
xt = 1-50 (dotted) and 1.34 (dashed) which are the 95% and 
99% goodness-of-fit of the models to the SN la data. 

we use the 54 data points from 'fit O of l]. For the CMB 
data the 99% goodness-of-fit limit is xt = 1-60 while the 
equivalent limit when using the SN la data is x 2 = 1-50. 
The figure clearly shows that smaller values of 7 are fa- 
vored and we can extract an upper limit to 7 of 7 < 1.6 at 
99% from the SN la data, 7 < 2.2 at 99% from the CMB 
data. This limit is a somewhat stringent because we have 
not marginalized over all the remaining parameters but 
the result is not very much different from other studies 
of IPL Quintessence such as Malquarti & Liddle 0] who 
derived 7 < 2 at 99%. Figure JHJ clearly demonstrates 
that if we use tracker initial conditions then the field must 
be very similar to the cosmological constant. Although 
we make this statement after considering only the IPL 
model this result is in qualitative agreement with other 
analysis of quintessence using the CMB and SN la which 
used a constant equation of state wq 0, |3^, 0, l4ll |42| . 



CMB, SN la without the tracker 

As we have seen from figure l |7| l. olt > 1 widens the 
allowed parameter space in 7 when fitting wq to present 
observed bounds. Larger values of ai behave more like 
a cosmological constant, wq = — 1, except for the caveat 
of an initial period when wq = +1. Thus, since larger 
values of ai push the current value of wq — > —1, we 
might expect that these solutions would provide a very 
good fit to the CMB and SN la data. 



10 



In figure © we plot the contours of constant \1 = 1.34 
and 1.50, that correspond to 95% and 99% goodness-of-fit 
to the SN la, and xl = 1-40 and 1.60, the 95% and 99% 
goodncss-of-fit to the CMB data, again in the j—a plane. 
For ai < 1 the exponent 7 is very tightly bound by the 
CMB and SN la because in this region the field joined 
the tracker solution before the present epoch. The limits 
on 7 are exactly the same as found for tracker initial 
conditions derived from figure ijSJ). But when 0.1 > 1 the 
field has not joined the tracker and we find the range of 
acceptable values of 7 increases significantly such that 
much of the region oli > 1 is compatible with SN la 
and CMB constraints. Thus the frequently stated limits 
on the exponent of power-law potential only apply when 
the initial value of ai falls within the range that will 
converge to the tracker solution by the present day for a 
given value of 7. 

V. CONCLUSIONS 

We have examined the IPL quintessence model focus- 
ing on the influence of the initial conditions upon the 
behavior of the solutions and in particular those cases 
where the field has not converged to the tracker solution 
by the present day. We have re-analyzed the previously 



derived bounds on the exponent 7 in the potential and 
find initial conditions that lead to solutions where the 
field reached the tracker before the present day limits the 
exponent to 7 < 2. However, after we looked at initial 
conditions where the field does not reach the tracker we 
found much larger values of 7 are permissible. The divid- 
ing line between the two types of solution occurs at the 
point of equipartition between the initial Quintessence 
and radiation energy densities and this is independent of 
the initial scale factor. IPL models of Quintessence give 
rise to values of the equations of state at the present time 
that are close to —1 either by the exponent being small, 
7 < 2, with small initial energy densities relative to the 
radiation or very large initial energy densities relative to 
the radiation content and arbitrary exponents. 
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